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ABSTRACT 

An  alternate  constitutive  formulation  for  visco-elastic  materials,  with 
particular  emphasis  on  macromolecular  viscoelastic  fluids,  is  presented  by 
generalizing  Maxwell's  idealized  separation  of  elastic  and  relaxation 
mechanisms.  The  notion  of  relative  rate  of  change  of  elastic  stress  is 
identified,  abstracted,  and  formulated  with  the  help  of  the  established  theory 
of  finitely  elastic  isotropic  materials.  This  gives  a  local  rate  type 
constitutive  relation  for  an  elastic  mechanism  in  a  simple  material. 

For  the  simplest  class  of  viscoelastic  polymer  melts,  the  notion  of  rate 
of  change  of  elastic  stress  and  its  damped  accumulation  is  identified  and 
formulated.  Under  conditions  of  moderate  strain  rates,  this  scheme  implies 
the  reliable  K-BKZ  model  for  a  class  of  polymer  melts.  An  obvious  extension 


generalizes  the  remaining  classical  spring-dashpot  models. 
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Nomenclature 


(1)  ^3  s  The  set  of  second  order  tensors  .  A-€tl  is  identified 

with  a  3x3  matrix  in  a  Cartesian  co-ordinate  system. 

(2)  0Sym:  The  set  of  symmetric  second  order  tensors. 


(3)  Q 


A  second  order  tensor  Q  is  orthogonal  if  and  only  if 
QT  -  Q-l. 
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This  is  a  symbol  for  the  functional 

SX  %3  sym»  where  X  is  the  set  of  piecewise 
continuous  and  differentiable  strain  histories 
JtQ:  U0,t]->X3  •  Other  functionals  appearing 
in°this  paper,  unless  otherwise  specified, 
should  be  interpreted  in  a  similar  manner. 
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ON  K-BKZ  AND  OTHER  VISCOELASTIC  MODELS  AS  CONTINUUM 
GENERALIZATIONS  OF  THE  CLASSICAL  SPRING-DASHPOT  MODELS 

* 

Amitabh  Narain 

1.  Introduction ;  The  K-BKZ  model  ([2],  [18],  [33])  for 
viscoelastic  fluids  is  known  to  give  good  experimental  fit  for 
rheometric  data  ([20],  [331)  obtained  for  an  identifiable  class  of 
melts.  However,  there  are  viscoelastic  fluids  and  strain  rates 
([21],  [38])  for  which  the  K-BKZ  model  does  not  give  good  results. 
Among  all  the  models  for  viscoelastic  fluids,  Maxwell's  model 
([28])  of  a  spring  and  a  dashpot  in  series  has  always  been  used  to 
give  a  qualitative  meaning  to  elasticity  and  relaxation  found  in 
these  fluids. 

In  this  paper  we  show  that  if  Maxwell's  notions  are  quantified 
and  generalized  for  an  idealized  strictly  viscoelastic  fluid 
possessing  only  the  two  mechanisms  of  fluid  elasticity  and  con¬ 
figurational  stress  relaxation,  then,  at  moderate  strain  rates,  it 
leads  to  a  K-BKZ  model.  The  single  integral  models  due  to 
Doi -Edwards  [13],  Curtiss  and  Bird  [12],  Bernstein,  Kearsley  and 
Zapas  [2],  Rivlin  and  Sawyers  [30]  are  closely  related  and,  in  a 
sense,  more  general  than  a  generalized  Maxwell  model.  A  survey 
article  of  Tanner  [333  and  rheometric  models  of  Wagner  ([36],  [37]) 
and  Laun  [20]  make  it  worthwhile  to  understand  the  scopes  and  limits 
tions  of  the  K-BKZ  model. 

Most  viscoelastic  models  fall  in  the  category  of  Noll's  simple 
materials  [25]  and  Coleman  and  Noll's  simple  fluids  [9].  Here  we 
abstract  and  generalize  the  essential  physical  notion  of  Maxwell  in 
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such  a  way  that  the  formulation  is  consistent  with  the  statements 
of  a  Simple  Material  (125),  [34]),  Objectivity  ([25],  [35]),  and 
material  symmetry  ([25],  [35])  associated  with  the  natural  rest 
states  of  a  simple  material. 

The  standard  definition  of  a  simple  fluid  is  given  by  Coleman 
and  Noll  [9]  in  terms  of  a  functional^  such  that  the  extra-stress 
is  given  by 

00 

T  +  pi  =  [Ct(t-s) ] ;  (1.1) 

£*0 

where  T  is  the  Cauchy  stress,  p  is  the  pressure,  and  Ct  (t-s)  is 
the  right  Cauchy-Green  tensor  evaluated  at  the  past  time  x  =  t-s<t. 

In  order  to  model  the  simplest  viscoelastic  fluid  within  the 
framework  of  (1.1),  the  notion  of  relative  rate  of  change  of 
elastic  stress  is  identified  in  the  established  theory  of  finitely 
elastic  isotropic  materials.  This  allows  changes  in  stress  on  a 
deforming  elastic  element  to  be  interpreted  as  a  response  to  the 
non-dissipative  changes  occurring  in  the  material  alignment  within 
the  element.  It  is  then  noted  that  a  change  in  material  alignment 
can  be  important  in  a  solid  as  well  as  in  some  fluids.  However,  a 
solid  has  a  preferred  material  alignment  in  its  rest  state  and  a 
fluid  has  no  preferred  material  alignment  in  its  rest  state.  This 
fact  is  used  in  introducing  and  formulating  the  concept  of  relative 
rate  of  change  of  elastic  stress  ^t(T)  for  an  incompressible  strict! 
viscoelastic  fluid.  A  strictly  viscoelastic  fluid  is  a  generalized 
Maxwell  fluid  endowed  with  the  two  mechanisms  of  elasticity  and 
configurational  relaxation.  Under  the  conditions  of  moderate 
strain  rates,  we  develop  the  following  generalized  Maxwell  model: 


T  =  -pi  +  /  G(°  )  ;t{t-o)do,  t>0 . 

!t<T)  =  It  i?(ct(T) > j f 

<MCt<T)>  =  1 4>i(  I ,  II )  Ct(T>  +  4>2(I,II)  Ct”l(T)]  (1.2) 

5  -k[  f ^( I ,11 )  Ct(T)  +  Ct  “ 1  <  t  )  J  , 


I  =  tr  Ct( t) , 

and 

II  =  tr  Ct_1  (x)  . 

In  (1.2),  G(0)  =  1,  Gf00)  =  0,  3  is  a  non-dimensional  monotoni- 
cally  decreasing  relaxation  function .  The  functions  and  f2  are 
scalar  valued  and  satisfy  the  conditions  of  zero  extra  stress  for  a 
rest  history 


fl  (3,3)  +  f2  (3,3)  *  0, 

and  the  condition  for  compatibility  with  linear  viscoelastic  fluids 


{ f  2 (  3  3)  -  f  ]_ (  3 , 3 ) )  + 


( 


•9fl 

3f, 

1 

jn 

(3,3) 

1 

91 

(3,3)] 

•8f2 

\ 

l3TT 

(3,3) 

z 

'  TT 

(  3 . 3  )J 

Once  a  Maxwell  model  is  generalized  to  (1.3),  we  present  obvious 
generalizations  of  other  spring  dashpot  models  for  Jeffrey's  liquid 
(see  [27])  and  the  standard  linear  viscoelastic  solid  ([14],  p.  175 
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2.  The  spring  and  dashpot  models. 


Me  will  now  recall  some  classical  spring  and  dashpot  models. 
For  a  long  time  it  has  been  recognized  (see  [11],  p.  26)  that,  if 
one  does  not  take  the  analytical  form  of  these  models  literally, 
the  physical  ideas  underlying  the  classical  mechanical  models 
should  have  a  valid  qualitative  generalization  to  general  three 
dimensional  straining. 

a.  A  spring  as  a  model  for  finite  elasticit 


When  a  non-linear  spring  (see  Fig.  2.1)  has  a  strain  e  (t)  with 
respect  to  its  unstretched  configuration  (t(to)=0)  occupied  at  some 


time  tQ,  the  stress  0(t)  is  given  by  a  function  g  such  that 
°(t)  *  g(£(t));  and  g(0)  *  0. 


(2.1 


V  k 
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Fig.  2.1:  A  spring  as  a  model  of  one  dimensional  elastic 
solid . 

This  constitutive  assumption  also  can  be  written  in  the  form 


ft  ■ 

O(t)  =  /  TT(X)dT, 


where 


5(1)  =  CT(r)  =  g(£  (t) ) ,T <  [tQ,t] . 


(2.2) 


This  emphasizes  the  notion  that  this  model  has  an  undamped  accumu- 
lation  of  elastic  stress  rate  n ( t ) .  We  wish  to  use  the  notion  of 


\V V-  ' 


accumulation  of  elastic  stress  rate  to  develop  a  constitutive 
theory  later  in  this  paper. 

b.  Maxwell  model  of  a  spring  and  a  dashpot  in  series . 


The  Maxwell  model,  with  a  spring  and  a  dashpot  in  series  ([27], 
p.  21),  was  meant  to  model  the  one  dimensional  straining  of  a 
viscoelastic  fluid  (see  Figure  2.2). 


Fig.  2.2:  A  schematic  representation  of  a  Maxwell  model. 

The  mechanical  devices  within  the  circle  represent 
physical  mechanisms  in  a  neighbourhood  of  the 
particle  P. 

If  e(t-s)  -  e(t)  is  the  relative  axial  strain  of  a  fluid  ele¬ 
ment  from  the  current  strain  e(t),  then  in  this  model,  the  stress  a 
on  the  particle  P  at  time  t  is  given  by 


—  00 

oo 

■L 


exp (-(t-t')/X)ke(t')dt' 


exp(-s/X)  {ke(t-s)}  ds , 


(2.3) 


where  e  is  the  strain,  and  X  =  ji.  In  Fig.  2.2,  k  is  the  linear 
spring  constant  and  n  is  the  dashpot  coefficient  ([14],  p.  175). 

Eq.  ( 2 . 3  >  2  has  the  following  alternate  interpretation. 
Principal  changes  in  stress  occur  through  the  elastic  mechanism 
represented  by  the  spring.  The  rate  of  change  of  elastic  stress 


Tt(t-s)  (or  at  any  past  time  x  =  t-s,  in  this  case,  is 


given  by 


Tt(t-s)  =  ke(t-s) 


(2.4) 


Because  of  the  relaxation  mechanism,  the  current  value  of 
stress  can  be  interpreted  as  a  damped  and  weighted  accumulation  of 
previous  elastic  changes  in  stress.  The  non-dimensional  relaxation 
function  exp(-s/A)  can  be  generalized  to  G(s)  and  a  Boltzmann  form 
of  the  Maxwell  model  may  be  written  in  the  form 


°(x,t)  =  J 


G(s)  Mt-s)  ds . 


(2.5) 


j 


In  (2.5),  G(s)  is  assumed  to  be  a  positive^  monotonically 
decreasing,  and  piecewise  continuously  differentiable  function  such 


that 


G(  0 )  =  1, 


(2.6) 


t*;. 


lim  G( s )  =  0 . 


The  above  interpretation  of  the  Maxwell  model  will  provide 
some  very  interesting  consequences  when  it  is  suitably  generalized 
later  in  Sections  5,  6,  and  7  of  this  paper, 
c .  The  standard  linear  viscoelastic  solid. 

The  standard  linear  viscoelastic  solid  ([14],  p.  175)  consists 
of  a  Maxwell  element  in  parallel  with  a  spring  (see  Pig.  2.3). 


•  ^  »  •  .  »  , 
•  *  '  * 

.V->V 


A  A 


Fig  2.3:  The  spring  in  the  above  model  provides  the 

solid  type  behaviour.  Viscoelastic  mechanism 
comes  through  the  Maxwell  component  of  the 
model . 

In  this  case,  it  is  well  known  that  a  in  Fig.  2.3  is  given  by 


o(P,t)  =  E£(t)  + 


exp(-s/>)  {ke(t-s)}  ds . 


(2.7) 


d .  Jeffrey's  model  with  Newtonian  viscosity. 

To  allow  for  noticeably  distinct  effects  of  rap  :.d  relaxation  at 
very  early  times,  presumably  associated  with  Newtonian  liquids,  it 
is  customary  to  allow  for  analogues  of  Jeffrey's  liquid  (27].  The 
mechanical  model  corresponds  to  a  Maxwell  element  and  ?  dashpot  in 
parallel  (see  Fig.  2.4). 


Fig.  2.4:  The  two  dashpots  have  the  coefficients  of 
and  r  .  The  linear  spring  constant  is  k. 

Here  it  follows  that  '  in  Fig.  2.4  is  given  by 

JO  •» 

'(P,t)  =  UTT  (  £  (  T  )  -  £  (  t)  )|  +  /  exp  (~s/‘)  { k  £  (  t  -5  )  }  ds  . 


3 .  A  brief  review  of  the  theory  of  simple  materials. 

Before  we  give  a  qualitative  continuum  generalization  of  the 
notions  in  section  2,  it  is  convenient  to  briefly  review  the 
general  scheme  of  Noll  [25]  and  Coleman  and  Noll  [9].  A  standard 
single  reference  for  such  matters  is  Truesdell  and  Noll  [35]. 

Let  X  be  the  place  occupied  by  a  particle  P  of  a  simple 
material  in  its  natural  stress  free  configuration  <  at  time  tQ. 

Let  the  current  position  of  the  particle  P  at  time  t  be  x,  and  let 
the  motion  for  T<[t0,t]  be  given  by  either  S  =  Xt  (x,?)  or  S  =Xt 
(X,f).  Clearly,  at  f  =  tQ,  we  must  have 

X  =  Xt(x,t0)  =  X tQ ( X  ,  t0 )  ,  (3.1) 

which  may  be  inverted  to  obtain 

x  =  Xt-1 ( X , t0 ) ,  (3.2) 

and ,  thus 

5-  Xt( Xt*1(X,t0) , T)  =  Xto(X,T). 

The  deformation  gradients  Fto  (X,t)  and  Ft  (x,t0)  are  defined  by 

Ft  (X,t)  =  Grad  Xt  (X,t), 

~'-o  -  '-o  - 

and  (3.3) 


Ft(x,t0)  =  grad  Xt(x,t0), 

and  we  note  that 

F  t  (x,tQ)  =  F  t  (X,t)_1.  (3.4) 

The  Cauchy  stress  T  (x,t)  of  a  homogeneous  simple  material  for 

which  the  deformation  is  measured  from  its  natural  configuration  • 

VC  w 

is  specified  by  a  response  functional  tv  ,  sucn  that 


Since 


T  =  M  <Ft  O). 


(3.5) 


8 


(3.6) 


ft  <2'T>  ftQ  <?*t)  *  ftQ  <*'T>  » 
it  then  follows  that 


K 

T  -  -H  (Ft<T>  FtQ  {t) >  •  (3.7) 

ic  *  t 

Note  that  the  form  of  the  functional  ]r\^  depends  on  the  natural 
configuration  <.  For  simple  materials  the  restriction  of 
Objectivity  ([25],  [34])  on  ,  is  given  by 


vc  ' 

Mv<Q(T)Ft(x)QT(t)  Q(t)Ft  (t)) 


=  Q(t) 


M, 


(Ft(T )Ft  (t) )QT(t) 


(3.8) 


for  any  orthogonal  tensor  function  Q(t),  tQ  <  t  <  t.  Using  the 
polar  decomposition  [16] 

Ft  (t)  =  Rt(r)  Ut(T)  -  Vt(T)  Rt(T)  (3.9) 

and  choosing 

Q(i)  »  Rt(x)T,  (3.10) 

we  then  see  that  (3.8)  yields 

K.  k 

H,(Ft(T)Fto(t))  =  H,  (0t<T)  Fto(t)).  (3.11) 


Thus ,  we  have 


T  =  KM,(Ut(T)  Ft  (t))  =  V(  ( U  t  (  t  )  ,  Ft  (t)). 


(3.12) 


The  new  response  functional  above  is  defined  and  introduced 

for  convenience.  Objectivity  requires  that  *"V(.  satisfy 


K  >  _  K. 


Q ( t )  H  <ut(t),  F  t  (t))  QT ( t )  =  Vi  ( Q ( t ) a  t ( T ) Qt ( t )  ,  Q ( t ) F  t  ( t ) ) 

'T.to  '  °  ~  ~  '  ~  ° 

( 3.13  ) 


for  every  orthogonal  Q.  The  form  of  (3.12)  emphasizes  that  only 
the  intermediate  local  shape  history  with  respect  to  the  current 


9 


configuration  (measured  through  Ut(T),  tD  <  t  <  t)  and  the  current 
local  shape  with  respect  to  the  reference  configuration  <  (measured 
through  FtQ(t))  affect  the  current  value  of  the  Cauchy  stress  at  a 
material  particle.  Now  the  particular  natural  rest  configuration  k 
in  the  development  above  may  be  indistinguishable  from  others  in  a 
set  of  natural  configurations  which  are  related  to  <  by  transfor¬ 
mations  belonging  to  a  symmetry  group  [16]  associated  with  the 
local  material  alignment .  A  common  statement  of  this  symmetry 
invariance  requirement  states  that  (c.f.,  [34]), 


T  =  H  <Ut<T>  >  Ft  (t)  H) 
~  ~  ° 

K  * 

=  M  (Ut(T)  ,  F*  (t> ) 
't-.Vo  ~  ~  ° 


(3.14) 


for  every  tensor  H  in  the  symmetry  group  /S  K  associated  with  k 
for  a  given  material.  For  an  isotropic  solid,  & K  is  the  full 
orthogonal  group  ([34]),  and  it  is  well  known  that  in  this  case 
(3.9)  and  (3.14)  may  be  combined  (with  H  =  RtQ(t)T)  to  give 


T  =  H  (UtU),  Vt  (t>>  . 
~  '  ° 


(3.15) 


One  finds,  using  the  uniqueness  of  the  decomposition  in  (3.9), 


that 


YtJt)-1  =  yt<tQ)  • 


(3.16) 


It  is  also  known  that  the  right  Cauchy  -  Green  tensor  Ct(x)  is 


given  by 


Ct(T)  =  ut(T)  Ut(  r)  ,t«  [tQ,t]  . 


(3.17  ) 


Ising  (3.16)  and  (3.17)  in  (3.15),  we  introduce  a  functional 


such  that 


.  -.1 

"•fi 

'•■'v'N'! 


'■  .  •’ 
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K.  ii 

T  =  ^  !Ct(T) ,  Ct(tQ) ] .  (3.18) 

^4,  ' 

In  (3.18),  the  material  is  in  its  undistorted  configuration  ug  to 
time  tQ  (i.e.,  C^It)  -  0  for  t  <  tQ) .  This  means  that  the  represent¬ 
ation  of  the  functional  in  (3.18)  does  not  depend  on  the  specific 
choice  of  the  undistorted  configuration  <  and  the  associated  time 
tQ  (see  Appendix  A).  This  allows  us  to  assert  that  the  natural 
state  of  an  isotropic  simple  solid  determines  a  response  functional 

\wsuch  that  i 

T  =  ^  rrvN^t(T)}.  (3.19) 

In  (3.19),  the  explicit  dependence  of  >3c^(4*ion  ct(to)  accounts  for 
the  differences  in  material  alignment  of  the  current  state  and  the 
preferred  isotropic  material  alignment  of  the  rest  state. 
Furthermore,  the  Objectivity  requirement  (3.13)  reduces  to 


■y  act<t0)Qr  toctiTitf1!  >Q  i  -y  &ia(?t«T))]?T 


(3.20) 


for  every  orthogonal  tensor  Q. 

For  an  incompressible  ([341,  p.  44)  simple  fluid, is  the 
full  unimodular  group  (i.e.,  any  ^  satisfies  det  H  =  1)  and  the 

resulting  constitutive  theory  reduces  to  (see  [34],  p.  57) 

t 

T  *  -pi  ♦  [Ct<  t)  J  ,  (3.21) 

with  the  constraint 

det  Ct(  t)  =  1,  T*[t0,t  J  . 

In  (3.21),  p  is  a  constituti vely  indeterminate  pressure  and  is 

independent  of  the  choice  of  reference  configuration.  Also,  (3.21) 
allows  tne  possibility  of  finding  an  undistorted  state  only  as  the 
time  tQ — >  -  1  for  a  given  0  <  t  <  ^.  However,  if  tQ  is  finite 


**  •  I 


.'•V  *1 

v.  >  .  • 
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(t0  >  -  “) ,  then  the  statement  (A. 12)  in  Appendix  A  allows  us  to 
set  tQ  ■  -  00  without  changing  the  value  of  the  stress.  Therefore 
without  loss  of  generality#  we  can  set  tQ  =  -  *>  in  (3.21)  and  thi 
makes  (3.21)  equivalent  to  the  statement  in  (1.1). 


4.  Finite  elasticity  and  the  notion  of  relative  rate  of  chanqe  of 


elastic  stress  St(T>. 

Let  us  consider  an  isotropic  solid  which  is  purely  elastic  and 
therefore  is  an  analogue  of  the  spring  in  section  2a.  In  this  case 
(3.18)  is  given  by 


T ( x , t )  =  K  ^3  [Ct(x),  Ct(t0)] 
=  fK  (Ct( tQ)) 


(4.1) 


where 


def  tc 

=  g  (B(t)), 


def  _ 

B(t)  =  FtQ(t)  Fto<t)T 

-  Ct(to)-1. 


(4.2) 


Q  g<  (B ( t ) )  QT  *  g<  (QB(t)QT) 

for  every  orthogonal  Q.  Similarly,  the  intermediate  stress  T(§,x) 
at  a  time  T-*[t0,t],  is  given  by  (see  Appendix  B): 

T ( 5 , t )  =  f <(CT (tQ) ) ,  (4*  3) 

*  $Ct(t0)  (Ct(T)), 

where  the  response  function  $Ct(t0)  depends  explicitly  on  Ct(tD). 
The  restriction  of  Objectivity  (see  (B.5)  in  Appendix  B)  becomes 


*QCt(t0)QT(Q  Ct(T)QT)  =  Q  [  $Ct(tQ) T) ) ^QT 


(4.4) 


for  every  orthogonal  tensor  Q  and  T4[t0,t], 

We  define  the  material  derivative  of  (4.3)  at  time  T  as  rate 


of  change  of  elastic  stress 


(  i)  and  write 


L*1  .* 1  f.  v.  i.",'  * 


T 


(C,T) 


d_ 

dx 


t*Ct(t0>  CCtCT>)l 


(4.5) 


*  nt  (x) . 

Since  at  time  tQ/  the  particle  is  at  X  and  the  material  element  is 
in  its  natural  configuration  <,  we  have 


T  (X,  tQ)  =  0. 

Collecting  (4.5),  (4.6),  and  ( B .  3 ) 3  from  Appendix  B,  we  have 


(4.6) 


r* . 

T  ( x ,  t )  =  /  TTt(x)dx 


-/ 


dT  l*Ct(t0)(Ct(T))]dx 


(4.7) 


=  $Ct(t0)<i>  “  9*  <B(t)). 

The  first  equality  in  (4.7)i  allows  us  to  interpret  the  current 
value  of  Cauchy  stress  as  an  undamped  accumulation  of  elastic 
stress  rate  This  stress  rate  measures  the  rate  of  non- 

dissipative  changes  occurring  in  the  material  alignment  of  the  ele¬ 
ment.  We  wish  to  borrow  this  notion  of  instantaneous  elasticity 
or  relative  rate  of  change  of  elastic  stress  TT^fx)  for  modeling 
instantaneous  elasticity  in  fluids.  Of  course,  for  a  fluid,  if 


tt*.  (  x ) 


-  S-[« 

dx 


ct(tQ)  (?t( 


x )  ^,  x 


-€  [tQ,tl 


(4.8) 


then  the  additional  restriction  of  fluid  symmetry  leading  to  (3.20) 
suggests  that  ^ct(tQ)  in  (4.8)  will  become  independent  of  Ct,(t0). 
This  is  because  the  rest  state  of  a  fluid  has  no  preferred  material 
alignment  but  there  are  differences  in  material  alignment  between 
two  unrelaxed  stressed  states  of  an  elastic  fluid  element. 


It  is  further  noted  that  it  is  impossible  to  have  an 
incompressible  purely  elastic  fluid  in  the  sense  of  an  elastic 


a 


solid  in  (4.7).  This  is  because  an  incompressible  purely  elastic 


fluid  will  have 


T  ( x ,  t ) 


/> 


(  T)dT 


=  $(1)  -  4>(Ct(t0) ) , 


(4.9) 


for  rrt(x  )  -  [  £  ( Ct  ( T  )  )  ]  . 

Now,  for  a  constant  Ct(tQ)  4  1  and  t  ^  T(x,t)  in  (4.9)  does  not 
tend  to  a  hydrostatic  pressure  for  an  arbitrary  isotropic  $ .  This 
is  in  contradiction  to  the  notion  of  a  rest  state  for  a  fluid. 

However,  as  we  see  in  the  next  sections,  a  fluid  (liJce  Maxwell 
model)  can  possess  instantaneous  elasticity  given  by  (4.9>2  if  it 
has  some  additional  mechanism  of  stress  relaxation. 


V-A 
V  V  V 
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5 .  On  phenomenological  motivation  for  modeling  a  class  of  polymer 


melts  as  a  generalized  Maxwell  model. 

We  now  wish  to  motivate  the  suggested  notion  of  rate  of  change 
of  elastic  stress.  For  this  purpose,  we  consider  a  class  of 
polymer  melts  as  a  representative  strictly  viscoelastic  fluid. 
Therefore,  in  what  follows,  the  following  specific  phenomenological 
picture  of  macromolecular  deformation  is  essential  to  the  validity 
of  the  modeling  presented  in  this  paper. 

Here  we  are  concerned  with  polymer  melts  which  contain 
numerous  supple  (flexible)  tubelike  strands  of  macromolecules  and 
their  networks.  During  deformation,  an  obvious  example  being 
extensional  flows,  an  arbitrary  material  element  of  the  fluid  may 
contain  somewhat  oriented  strands  of  macromolecules.  Yet  the 
stress  free  prolonged  rest  state  corresponds  to  tensionless  strands 
of  macromolecules  which  are  randomly  oriented  within  a  specified 
fluid  element.  This  random  orientation  of  inert  tensionless  macro¬ 
molecules  is  required  for  consistency  with  the  unimodular  symmetry 
[34]  expected  from  the  prolonged  rest  state  of  any  simple  fluid. 
However,  when  a  stressed  fluid  element  is  deformed  from  one  shape 
to  another,  the  supple  resilient  macromolecules  conform  to  the  new 
shape  and  some  stress  is  generated  through  tension  in  the  macromo¬ 
lecules.  As  the  fluid  element  is  made  to  change  its  shape  by 
incremental  stress,  this  new  deformed  shape  has  some  elastic  pro¬ 
pensity  to  return  to  its  original  shape.  Accompanying  this  moment¬ 
ary  non-dissipative  change  in  the  alignment  of  the  macromolecules, 
on  a  lagging  time  scale,  there  is  a  release  of  tension  through 
dientanglement  and  readjustment  of  macromolecular  networks.  This 


release  of  tension  within  an  oriented  shape  of  macromolecules  is 
assumed  here  to  be  the  principle  mechanism  of  configurational 
relaxation .  This  notion  is  perhaps  the  same  as  the  notion  of 
configurational  relaxation  proposed  by  Giesekus  [15].  Of  course, 
if  the  fluid  element  stops  deforming,  this  configurational  relaxa¬ 
tion  will  eventually  produce  a  state  of  tensionlessness  and  this 
state  will  easily  allow  randomizing  of  macromolecular  orientation 
to  take  place.  Therefore,  any  randomizing  associated  with  relaxa¬ 
tion  is  assumed  to  be  noticeable  only  in  cessational  flows 
(straining  of  a  fluid  element  stops),  and  even  in  this  case  config u 
rational  relaxation  is  assumed  to  be  the  dominant  part  of  relaxa¬ 
tion.  Therefore,  we  assume  that  changes  in  stress  are  due  to  two 
mechanisms:  one  is  through  an  elastic  mechanism  and  the  other  is 
through  a  relaxation  mechanism.  The  elastic  mechanism  corresponds 
to  resilient  alignment  of  supple  macromolecules  within  the  changing 
shape  of  a  fluid  element.  We  further  assume,  despite  the  presence 
of  some  possible  reptations  [13],  that  all  effects  of  material 
alignment  can  b^modeled  by  the  changing  shape  of  a  fluid  element. 
Since  we  are  interested  in  the  current  value  of  the  Cauchy  stress  . 
it  is  beneficial  to  use  the  current  configuration  of  a  fluid  ele¬ 
ment  for  modeling  different  elastic  contributions  from  differently 
aligned  macromolecules  in  shapes  different  than  the  current  one. 

For  this  class  of  polymer  melts,  as  said  earlier,  we  assume  that 
the  effects  of  the  tendency  of  macromolecules  to  randomize  within  a 
deforming  material  element  is  insignificant  when  compared  to 
configurational  relaxation  and  the  elastic  mechanism  associated 
with  direct  mechanical  straining*.  This  does  not  mean  that  all 

*The  fact  that  randomization  in  some  melts  and  concentrated  solu- 
tions  can  be  ignored  is  also  discussed  by  Marrucci  (Polymers, 
Liquid  Crystals,  and  Low  Dimensional  Solids,  Plenum  Press,  d 
149 ,  1984  )  .  p 


polymers  can  be  modeled  by  ignoring  effects  of  macromolecular  ran¬ 
domization  within  a  given  shape  of  a  material  element.  In  fact, 
the  known  theories  ([5],  16J)  of  rubber  elasticity  and  other 
theories  of  entropic  elasticity  (111,  [6]),  can  possibly  be  adapted 
into  this  proposed  modeling  approach  to  include  the  effects  of 
macromolecular  randomization.  However j  in  this  paper,  we  restrict 
ourselves  to  those  melts  which  have  ignorable  contribution  asso¬ 
ciated  with  randomization. 

Therefore,  when  a  strictly  viscoelastic  element  (identified  at 

def 

the  current  time  t)  is  strained  from  a  configuration  C(t)(o+Aa)  * 

ct(t  -(a+Aa))  to  a  configuration  C<t)  (^),  the  change  in  shape  acti- 

Ao 

vates  a  direct  change  in  stress  (modeled  as  elastic  change)  * 

• 

Tt  (t-o)Acr+  o  (Aa)  near  the  past  time  t  *  t-c  .  Having  looked  at 
the  notion  of  relative  rate  of  change  of  elastic  stress  in  section 
4,  we  borrow  the  constitutive  model  (4.8)  for  dependence  of  "^(t-o) 
on  the  strain  Ct(t-o)  and  on  the  strain  rate  Ct  (t-o). 

Because  of  the  presence  of  relaxation  mechanism,  the  current 
value  of  Cauchy  stress  T(x,t)  is  a  damped  accumulation  of  incremen¬ 
tal  stresses  This  notion  is  a  qualitative  generalization  of 

Maxwell  model  and  is  formally  stated  in  the  next  section. 


On  idealized  presence  of  elastic  and  dissipative  mechanisms  for 
a  generalized  Maxwell  model. 

For  a  strictly  viscoelastic  fluid,  we  assume  the  standard 


kinematics  of  section  3.  In  addition,  at  any  instant  of  time  t, 
while  chasing  an  identifiable  neighbourhood  of  a  particle  P,  there 
exists  a  rate  of  change  of  elastic  stress  history 

^(t):  10,»)  -  ->  t3sym,  (6.1) 

•  JL  • 

and  71  ( )  (  s  )  =  MP,t-s)  =  rj.(x,t;t-s)  is  the  relative  rate  of 

change  of  elastic  stress.  For  x^[tQ,t]  and  tg  >  —  00  ,  (4.8)  is 
written  as 

!t<T)  =  It  t$Ct(t0)  <Ct(T ;) )  1 

(6.2) 

=  D*  Ct(tQ)  (Ct(T))  Ct(T). 

In  (6.2),  the  strain  dependent  elasticity  modulus  "P'l’Ct  ( tQ )  (Ct(x))" 

models  the  effect  of  material  alignment  in  changing  shapes  and  the 
strain  rate  Ct ( x )  is  the  principle  cause  of  changes  in  stress. 

For  now,  we  postulate  tt  t  ( x )  to  be  given  by  (6.2)  subject  to 
further  restrictions  imposed  by  the  increased  symmetry  of  the  rest 
state  of  a  fluid.  To  model  the  effect  of  relaxation  and  its  possible 
coupling  with  elasticity,  we  postulate  that  the  current  value  of  extra 
stress  for  an  incompressible  generalized  Maxwell  fluid  is  given  by 

T  (  x  ,  t )  =  -pi  +  I  TTt  (  x  )  J  ,  (6.3) 

~  'T.k,  ' 

where  is  an  accumulation  functional  associated  with  relaxation. 
Without  loss  of  generality,  the  rest  state  is  taken  to  be  at  t0  =  -  ■ 


and  (6.3)  is  rewritten  as 


T  -  -pi  +  *H  t  TTfc  (  t-S  )  ]  =■  -pi  +  l"(t>  1  • 

-  '  s«o 


(6.4) 


On  substituting  (6.2)  in  (6.4)  and  then  applying  the  notion  of 
fluid  symmetry  leading  to  (3.21),  one  can  show  that 


*t<T>  =  T7  [*(Ct(T) ) ] , 


(6.5) 


for  T4(-°°,t].  Clearly  (6.4)  and  (6.5)  together  give  a  special 
class  of  simple  fluid  given  by  (1.1).  The  mathematical  assumptions 
on which  leads  to  (6.4)  and  (6.5)  may  be  an  interesting  area  of 
investigation.  However,  in  this  paper,  as  is  the  case  with  the 
original  Maxwell  model,  the  notion  of  elasticity  and  relaxation  are 
modeled  directly  on  the  basis  of  the  phenomenology  described  and 
assumed  in  section  5. 

We  now  see  that  the  restriction  of  Objectivity  in  (4.4)  and  (4.8) 


implies  isotropy  of  the  function  <f>  in  (6.5).  Therefore 
4  (Q  Ct(T)  qT)  =  Q  4> ( Ct (  T )  )  QT 


(6.6) 


for  any  orthogonal  Q.  It  is  noted  that  the  restriction  (6.6)  is 
imposed  by  the  fact  that  our  Euclidean  space  has  no  preferred 


directions  and  it  has  nothing  to  do  with  material  isotropy  at  the 


strain  Ct(T)*  In  fact,  we  recall: 

C(Ct(t-s) )  *  D$  (Ct(t-s) ) 


(6.7) 


in  (6.2)  measures  the  effects  of  relative  differences  in  material 
alignments  between  the  continuum  configurations  at  times  t  =  t-s 
and  ?  =  t.  It  is  instructive  to  relate  the  choice  in  (6.5)  to  the 
notion  of  elastic  resilience  alluded  to  in  the  phenomenological 
description  of  section  5.  For  this,  consider  a  past  time  r  =  t-(.'+\ ') 


and  the  strain  Ct  (  :  -  (  '+•  \  • )  )  =  C<  t )  (  h.\o  )  .  When  the  fluid  element 


is  brought  to  the  configuration  C(t)(o),  we  do  so  by  an  incremental 


stress 


t-o 


Ao 

"t 


■/ 


TTt  (  t-S)dS 

t- (o+Ao) 

*  £  4> C Ct ( t-a) )  -  4>(Ct(t~(  a+Aa) )  ] 

~  <V  «v  «s» 

*  $(C(t)(a»-  *(C(t) < a+Aa) ) . 

^  •* 

Now  before  relaxation  sets  in,  if  the  incremental  stress 
set  to  zero,  (6.8)  requires  that  we  must  have 


(6.8) 


Ao 


was 


!(C(t)  (  0+Aa)  >  "  $(Q(t)(o))  =  0.  (6.9) 

In  the  thought  experiment  resulting  in  (6.9),  we  find  that  if  <*>  is 
invertible  at  C(t )  <  a)  (that  is  the  elastic  modulus  D<J>  (C(t)(o))  is 
such  that  D<J>  (C(t)(o))A  =  0  implies  A  =  sym)  ,  then  (6.9) 

implies 


C(t)(a+Ao)  =  C(t)(a).  (6.10) 

We  view  (6.10)  to  be  in  accord  with  the  physical  notion  of  elastic 
recoil  (see  [4])  when  the  stress  is  removed  prior  to  the  onset 

of  relaxation. 

sup 

Furthermore,  if  the  maximum  magnitude  of  strain  t  ( .ro(t  j 
I  ?t*T*“i  !  exists  and  is  close  to  zero,  then  ^(t-s)  in  (6.7)  can 
be  approximated  by 

5(t) <»)  =  ®  Ct ( t-s ) ,  (6.11) 

where  C  =  D<5  (1)  is  a  constant  fourth  order  tensor  in  the  space  of 
linear  maps  from  sym  to  ^symt  The  constant  tensor  C  in  (6.11) 
makes  sense  because  now  there  is  no  significant  difference  in 
material  alignments  at  all  times  r  <  t.  On  applying  the  isotropy 
condition  (6.2),  we  find  that 

Q  1C  Ct(t-s)]  Qt  =  C  [Q  Ct ( t-s )  QT]  (6.12) 

for  any  orthogonal  tensor  Q.  Using  a  well  known  result  [32)  often 


used  in  linear  theory  of  elasticity,  it  follows  that  (6.12)  implies 
that  C  Ct(t-s)  must  be  of  the  form 

X  Ct ( t-s )  =  k  Ct(t-s)  +  v  tr(Ct(t-s))  1  (6.13) 

for  some  constants  k  and  v.  Eq.  (6.12)  and  (6.13)  embody  the  notion 
of  a  linear  elastic  mechanism  and  correspond  to  a  linear  spring  in 
the  Maxwell  model  of  Fig.  2.2.  The  general  isotropy  condition  in 
(6.6)  along  with  an  adaptation  of  theroems  of  Rivlin  and  Ericksen 
[29],  and  Serrin  [32]  implies  that  $  has  a  representation 

®(Ct(t -a))  =  4>0 ( 1 ,11 )  1  +  <f>2.( 1 ,11 ) Ct ( t-o)  +  4>2<I  ,11 ) Ct-1  ( t-o )  , 
where  (6.14) 

I  =  tr ( Ct ( t-a ) ) ,  and  II  =  tr ( Ct_1 ( t-o ) ) . 

As  a  result  of  the  incompressibility  condition  (3.21)2#  the  third 
invariant  of  Cc(t-o)  does  not  appear  in  (6.14).  <J>q,  $i#  and  4> 2 
eq .  (6.14)  are  material  functions  of  the  invariants  I  and  II. 
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tion  h  be  small.  Of  course  we  assume  that  the  norm  of  0$  remains 
bounded  in  the  space  of  linear  maps  from  ^sYm  to  ^JsYm.  Therefore, 
this  assumption  allows  for  large  and  even  unbounded  magnitudes  of 
strains  as  long  as  Ct(t-s)  is  small  in  the  recent  past.  We  now 
make  an  additional  physical  assumption  regarding  G.  Let  G  be  posi¬ 
tive  on  [O,00),  piecewise  continuously  differentiable  on  (0.°°), 
monotonically  decreasing  on  lO,00), 

G  ( 0 )  =  1, 

and  (7.5) 

G  (  s  )  — ►  0  as  s  — *  oo 

sufficiently  fast  for  a  meaningful  integral  in  (7.3).  The  assump¬ 
tion  (7.5)  on  G  in  (7.3)  can  be  shown  to  follow  from  requiring  that 
the  magnitude  of  the  traction  vector  on  the  boundary  of  a  fluid 
element  must  never  increase  if  the  element  was  made  to  experience  a 
static  continuation  (see  [26]  for  the  meaning  of  this  term). 

We  note  that  ^(t)  has  the  physical  dimensions  (ML'^-T-3)  of 
rate  of  change  of  stress.  This  implies  that  G(s)  is  nondimen- 
sional.  This,  in  turn,  suggests  a  provable  implication  that  there 
exists  at  least  one  relaxation  time  T  and  a  function  G  such  that 

6(s)  •  G(s/X)  (7.6) 

for  every  s_>  0. 

If  we  now  consider  histories  C(^)  which  are  such  that  both 
I  I  C(t)-1  |  |  ft  and  |  |  C(t)  |  |  h  are  small  in  the  inner  product  space 

induced  by  (7.1),  then,  the  linearization  (7.3)  should  be  expected 
to  reduce  to  the  standard  and  general  result  of  Coleman  and  Noll  in 
(5.18)  of  [10].  For  such  a  history  we  use  the  linear  relation  for 
" ( t )  in  (6.13)  and  substitute  it  in  (7.3). 
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After  absorbing  the  multiples  of  unit  tensor  in  the  indeterminate 
pressure,  we  find 


-jLin  =  -pi  +  /  S(s)  k  Ct(t-s)ds 


(7.7) 


The  result  in  (7.7)  is  compatible  with  the  standard  linearization 
(cf.  (5.18)  of  110]) 


ijLin  s  -pi  +J  G'(s)  Ct,(t-s)ds. 


(7.8) 


The  function  G  and  G  in  (7.7)  and  (7.8)  are  related  by 
G  ( 0 )  =  k  >  0  , 


def 

G(s)  =  k  G(s) ,  s>  0 . 


(7.9) 


The  result  in  (7.9)  gives  an  important  indication  that  the  postu¬ 
late  in  (6.3)  is  a  good  one.  This  is  because  ^G(0)/p  is  known  to 
give  the  shock  speed  ([8],  [23))  for  a  linearized  strictly 
viscoelastic  fluid  given  by  (7.8).  Interpretation  of  G(0)=k  as  ' 
modulus  of  the  elastic  mechanism  in  the  fluid  seems  to  give  meaning 
to  this  result.  The  shock  speed  is  therefore  due  to  elastic  resi¬ 
lience  to  impact.  Also  various  results  of  Coleman,  Gurtin  and 
Herrera  [9],  Narain  and  Joseph  [23],  and  Renardy  [28]  show  that  the 
amplitude  of  a  shock  decays  with  time  as  exp  ( tG' ( 0 )/G( 0 ) ) .  In  the 
above  interpretation,  the  decay  is  governed  only  by  the  parameter 
S'(0)  =  G'(0)/G(0)  <  0,  and  this  parameter  is  related  only  to  the 
relaxation  mechanism  of  a  linearized  strictly  viscoelastic  fluid. 
Similar  interpretation  of  other  formulas  of  Narain  and  Joseph  [23] 
is  discussed  in  Narain  [22]. 


8.  The  assumptions  underlying  the  K-BKZ  type  models  for  polymer 
melts . 

In  a  survey  article/  Tanner  [33J  finds  that  the  K-BKZ  type 
model  is  superior  regarding  compatibility  with  experiments  in  a 
large  class  of  motions  of  polymeric  fluids.  Tanner  [33]  indicated 
that  there  is  no  good  qualitative  understanding  for  this  model.  It 
is  hoped  that  this  paper  contributes  towards  this  purpose.  The 
representation  of  the  K-BKZ  model  seems  to  be  based  on  the  assump¬ 
tion  that  the  fluid  can  be  modeled  as  a  generalized  Maxwell  model 
and  that  the  strain-rates  are  moderate  enough  to  ignore  the  second 
term  in  (7 . 3  >  i .  On  substituting  (6.5)  in  (7 . 3 )  r  and  ignoring  the 
smaller  term,  one  finds 


T 


d_ 

ds 


$(Ct (t-s) )ds . 


(8.1) 


The  assumption  of  moderate  strain  rates  seems  to  be  justified  by 
the  fact  that  the  K-BKZ  model  (8.1)  is  not  accurate  (see  [21]) 
under  high  strain  rates  involved  in  experiments  with  multiple  steps 
in  strain  .  This  is  because,  under  these  conditions,  |  |  Ct  j  |  h  and 
hence  o  (  |  |  I  !  h*  *-n  (7.3)].  is  not  negligible.  Physically,  this 
amounts  to  saying  that  the  time  scale  associated  with  rate  of 

strain  is  small  compared  to  a  characteristic  time  scale  X  of 

-  U 

relaxation  (high  Weissenberg  Number,  X  j-^see  (3.15)  of  [11]).  That 
is,  rate  of  elastic  energy  build  up  far  exceeds  the  dissipation 
rate  and  therefore  the  assumption  of  independent  elastic  and 
relaxation  mechanisms  in  (8.1)  is  not  very  good. 

On  substituting  (6.14)  in  (8.1)  and  integrating  it  by  parts 


while  absorbing  the  multiples  or  unit  tensor  1  in  the  pressure 
term,  we  find 


T  =  -pi  •I  ?Mo)(^|I(II)Ct(t-o)  +  $2(I  ,II>Ct*1(t-o)}  da.  (8. 


2) 


One  can  relate  (8.2)  to  the  popular  K-BKZ  form  ([33]),  by  intro¬ 
ducing 


def 


a(°)  -  -k  G' (o)  »  -G' (O) fo  >  0 


1  def 

£  $1(1,1!)  *  -fi(I, II) , 


and 


(8.3) 


1  def 

£  $2(I'I;I)  *  -f2<I»II) 


where  k>0  is  the  same  as  in  (7.7). 

Using  (8.3),  (8.2)  is  written  in  the  form 


T  »  -pi  +  f  a(0){f1(i,ll)Ct(t-o)  +  f 2(I ,II)Ct_1(t-a)}d0. 


(8.4) 


The  functions  fi  and  f2  in  (8.4)  can  be  normalized  by  requiring  the 
extra-stress  to  be  zero  for  the  rest  history  Ct(x)  =  1,  r  <  t. 

This  gives 

f  1  (  3 , 3 )  +  f  2  (  3 , 3 )  =  0. 

Furthermore,  a  linearization  of  (8.4)  gives 


(8.5) 


T  =  -pi  +  /  a (cr)  C  Ct(t-o )do , 
0 


(8.6) 


where 


Pfl 

=\—  <3' 


3f, 


3)  - 
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(3 


l 

I3f2 

9f-  3 

>J 

{3'3> 

-  m  °'3,J 

+  { f  1  (  3 , 3  )  -  f  2  (  3  ,  3  )  }  . 


For  (8.6)  to  be  compatible  with  (7.8),  it  is  required  that 


8.  The  assumptions  underlying  the  K-BKZ  type  models  for  polymer 
melts . 

In  a  survey  article.  Tanner  [33]  finds  that  the  K-BKZ  type 
model  is  superior  regarding  compatibility  with  experiments  in  a 
large  class  of  motions  of  polymeric  fluids.  Tanner  [33]  indicated 
that  there  is  no  good  qualitative  understanding  for  this  model.  It 
is  hoped  that  this  paper  contributes  towards  this  purpose.  The 
representation  of  the  K-BKZ  model  seems  to  be  based  on  the  assump¬ 
tion  that  the  fluid  can  be  modeled  as  a  generalized  Maxwell  model 
and  that  the  strain-rates  are  moderate  enough  to  ignore  the  second 
term  in  (7.3)^.  On  substituting  (6.5)  in  (7.3>i,  and  ignoring  the 
smaller  term,  one  finds 


-pi  +  j  -G(s) 
JQ 


<MCt(t-s)  )ds. 


(8.1) 


The  assumption  of  moderate  strain  rates  seems  to  be  justified  by 

the  fact  that  the  K-BKZ  model  (8.1)  is  not  accurate  (see  [21]) 

under  high  strain  rates  involved  in  experiments  with  multiple  steps 

in  strain  .  This  is  because,  under  these  conditions,  j  j  Ct  |  |  h  anc* 

hence  o  (  |  I  ^t  I  I  h5  *n  <7*3>i  is  not  negligible.  Physically,  this 

amounts  to  saying  that  the  time  scale  associated  with  rate  of 

strain  is  small  compared  to  a  characteristic  time  scale  \  of 

-  U 

relaxation  (high  Weissenberg  Number,  X  £:see  (3.15)  of  [11]).  That 
is,  rate  of  elastic  energy  build  up  far  exceeds  the  dissipation 
rate  and  therefore  the  assumption  of  independent  elastic  and 
relaxation  mechanisms  in  (8.1)  is  not  very  good. 

On  substituting  (6.14)  in  (8.1)  and  integrating  it  by  parts 
while  absorbing  the  multiples  of  unit  tensor  1  in  the  pressure 
term,  we  find 


*  ^  • 

.  •  .  *  » 


*■  .  •  «  *  « 


*wfcv. 


BuL. ... 


ate 


21 


-  «  .  -  _  *  „  «r  .  ■  _  t  r„  « 


r 00 

T  a  -pi  +/  ?'(0){^(l,II)Ct(t-0)  +  02(I do-  (8* 

-  J  o 

One  can  relate  (8.2)  to  the  popular  K-BKZ  form  (1331),  by  intro¬ 
ducing 

def 

a (o )  -  -k  G' (a)  =  -G' (a) ,o  >  o 


2) 


and 


i  def 

£  ♦iiMi)  =  -fid  ,ii) , 


i  def 

i  0  2  ( 1 ,11)  =  -f2(I»H)  » 


(8.3) 


where  k>0  is  the  same  as  in  (7.7). 

Using  (8.3),  (8.2)  is  written  in  the  form 


T  =  -pi  +  J  a<a)  {f!(I  ,II)Ct(t-o)  +  f  2  (I  ,H  )  Ct_1  ( t-o )}  do  .  (8. 


4) 


The  functions  fi  and  f2  in  (8.4)  can  be  normalized  by  requiring  the 
extra-stress  to  be  zero  for  the  rest  history  Ct(t)  =  1,  t  <  t. 

This  gives 

f 1 ( 3 , 3 )  +  f  2 ( 3 , 3 )  =  0. 

Furthermore,  a  linearization  of  (8.4)  gives 


(8.5) 


T  ■  “Pi  +  /  a(o )  C  Ct<t-a )do , 
0 


(8.6) 


where 


"131 


3f- 


31  (3,3)  -  3ZI ( 3 


A  -ft 


3f,  3f, 

1  (3,3)  -  jjf  (3 


,3)} 


+  { f i ( 3 , 3 )  -  f 2(3,3)  }. 

For  (8.6)  to  be  compatiDle  with  (7.8),  it  is  required  that 

C  =  -i. 
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The  suitability  of  the  above  K-BKZ  type  model  can  only  be 
verified  by  applications  and  experiments.  Important  semi -empirical 
formulations  of  Wagner  ([36],  [37]),  Laun  [20],  and  Chang  and  Lodge 
[7]  can  easily  be  shown  to  lie  in  the  above  framework  with  f].(£,II)»0. 

Since  Lodge’s  rubberlike  fluid  [3]  is  a  special  case  of  (8.4), 
various  experimental  comparisons  [33]  suggest  that  it  is  advisable 
to  do  melt  rheometry  in  this  more  general  framework.  The  most 
important  experimental  verification  and  suggestion  of  the  model  in 
(8.4)  can  be  found  in  Laun ’ s  paper  [20].  This  is  briefly  rein¬ 
terpreted  in  section  11 . 


In  this  section  we  generalize  the  notion  of  a  standard  linear 
viscoelastic  solid  shown  in  Fig.  2.3. 

For  an  isotropic  viscoelastic  solid,  the  development  of  sec¬ 
tion  3  leading  to  (3.14)  holds.  Using  B(t)  in  (4. 1)4  and  the  fact 
that 

B(t)  =  vto  (t)  yto(t> 

uniquely  determines  VtQ(t)  for  a  given  B(t),  one  can  replace 
in  (3.15)  by  a  functional such  that 


T(x,t)  =  >(  (Ut(T )  ,  Vt  (t) ) 

~  ~  0 

=  (B ( t )  ,  Ct(T  )  )  . 


(9.1) 


Without  loss  of  generality,  we  can  set  t0  =  -  00  in  (9.1)  to 
rewrite  it  in  the  more  familiar  form 

K.  i- 

T  =  (B(t)  ,  Ct(T ) >  (9.2) 

'Ts-O*"' 

given  by  Truesdell  and  Noll  (35J.  We  wish  to  seek  a  specific 
approximate  representation  of  (9.2)  for  a  generalization,  under 
conditions  of  moderate  strain  rates, of  the  model  in  Fig.  2.3. 

For  this  we  imagine  a  homogeneous  matrix  of  a  viscoelastic 
solid  with  a  stress  free  natural  configuration  <  exhibiting  isotro¬ 
pic  material  symmetry.  It  is  helpful  to  imagine  that  this  matrix 
is  made  up  of  random  distribution  of  elastic  rod-like  solids  and 
r .  <■>  x  i  b  1  c  melt  like  macromolecules. 

During  any  straining  at  a  time  [ t0 ,  tj,  the  relative  rate 
change  of  elastic  stress  ” t ( x )  with  respect  to  the  current  con¬ 
figuration  at  time  t  can  be  modeled  as 


30 


where 


Ht < T)  =  Ut  < T)  +  *t( T) » 


it<T>  -  §?  t*ct<t0><ct(T>>] 


is  the  same  as  in  (4.5)  and 


?t<T>  =  If  I  *  c  Ct  C  T ) ) ] 


(9.3) 


(9.4) 


(9.5) 


is  the  same  as  in  (6.5).  Furthermore,  by  analogy  to  the  case  in 
Fig.  2.3,  we  can  postulate  that  the  evolution  of  the  Cauchy  stress 


at  any  T4tt0,t]  is  given  as 


where 


T  ( C, i)  ■  Ts( €,t>  +  Tf  ( 5,t) , 


Ts  (C,t)  =  TTt(C,T) 


(9.6) 


(9.7) 


T 

if  ‘i'T)  =  • 


Now  integration  of  (9.6)  for  an  incompressible  solid  yields 

V 

T ( x , t )  =  -pi  +  gK  ( B ( t ) )  +  t^t  (T)J- 


(9.8) 


(9.9) 


The  function  g  in  (9.9)  is  as  defined  in  (4.1).  Let  tQ  — >  -  03  and 
B ( t )  =  £1™_00  be  such  that  the  solid  part  of  the  matrix  is 

in  its  elastic  realm  governed  by  (9.7).  Following  the  approxima¬ 
tion  leading  to  (8.4),  at  moderate  strain  rates,  T  in  (9.9)  can  be 
written  as 


T  =  -pi  +  g  ( B ( t ) )  +  /  ( -G 1 ( a ) )  { f x ( I ( t-o ) , I I ( t-a ) )  Ct ( t-' )  + 


f2  ( I  ( t  -  : )  ,  1 1  (t-.-nCt'1  (t-o)f  do. 


(9.10) 


where 


g  ( B  ( t )  )  =  •- 1  <  1 1  ( t )  ,  1 2  ( t )  )  S(t) 


(9  .  il  ) 


The  last  term  in 
and  62  in  (9.11) 


+  h  dlHt)*  X  2  <  t- )  > 

(9.10)  is  as  in  (8.4). 

are  functions  of 
def 

ll<t)  ■  trB(t), 


B“i  ( t )  . 

The  material  functions  Bj_ 

(9.12) 


and 


def 


I^(t) 


trB"i(t) . 


10 .  The  generalization  of  Jeffrey's  liquid,  and  Kelvin-Voigt  Solid 
For  polymer  solutions  in  a  low  molecular  weight  Newtonian 
solvent,  one  expects  that  the  effect  of  the  solvent  is  a  rapid 
relaxation  associated  with  the  dashpot  in  Fig.  2.4.  The  first  term 
in  the  right  side  of  (2.8)  can  be  shown  to  generalize  to  a 
Newtonian  term 

UgY  (G(T)-e(t))|  Tat->uCt(I)  I  T=t  (10.1) 

*  2u  D(x ,t ) . 

In  (10.1)  above,  D  is  the  symmetric  part  of  tne  velocity  gradient 
"grad  V(x,t)“.  Using  (10.1)  and  (8.4)  it  follows  that  (2.8)  for  an 
incompressible  viscoelastic  fluid  may  be  generalized  as 

T  *  pl+2MD+  /  (-G' (a) ){f1(I,II)Ct(t-o)+f2(I,II)Ct-i(t-a)}da.  (10.2) 

^0 

It  is  well  known  that  (10.2)  also  follows  from  the  strictly  visco¬ 
elastic  model  in  (8.4)  if  we  formally  replace  G(s)  by 

u6 ( s-0 )  +  G(s) ,  (10.3) 

where  6 (s-0)  is  the  Oirac  delta  function  and  u  >0.  A  theoretical 
justification  for  allowing  distributions  in  the  kernel  of  the 
linearized  constitutive  relation  has  been  discussed  by  Saut  and 
Joseph  [31J.  It  is  conceivable  that  the  so  called  exact  Newtonian 
contribution  in  (10.3)  is  merely  an  idealized  way  of  asserting  a 
rapid  relaxation  present  in  a  regular  G(s)  near  s=0 . 

Accepting  the  generalization  in  (10.1)  and  (4.7),  a 
Kelvin-Voigt  incompressible  solid  ([17],  p.  21)  generalizes  to  give 
T  =  -pi  +  g<  (Bit))  +  2  D  ( x  ,  t )  .  (10.4) 

Similarly,  a  parallel  combination  of  a  Kelvin-Voigt  element 
and  a  Maxwell  element  generalizes  to  give 


11 .  Compatibility  of  the  K-BKZ  model  with  some  successful 
semi -empirical  models  in  shear  flows. 

Of  particular  significance  is  the  fact  that  some  successful 
semi -empirical  models  in  shearing  motion  are  implied  by  the  K-BKZ 
model.  Consider  a  simple  shearing  motion  in  y  direction  and  in 
x-y  plane;  the  motion  x  ■  x(X,t)  of  a  particle  X  (X,Y,Z)  is 

•V 

given  in  Cartesian  co-ordinates  by  the  relation 

(x,y,z)  *  (X,  Y  +  y(x,t),  Z).  (11.1) 

The  shear  strain  Y(t)  and  the  relative  shear  strain  yt(s)  are 
defined  as 


3ji 

Y(t)  » 


and 


3u  3  u 

Yt<8^  “  t-s)  -  (x,t) 


(11.2) 


(Ct(t-s) l*[C(t j (s) ] 


A  standard  calculation  shows  that  the  Cartesian  components  of 
C(t) (s)  is  given  by 

r  2 

1+Y  t(sP  Yt(s)  0 

Yt { s )  1  0 

0  0  1 

The  corresponding  components  of  Ct-1(t-s)  is  given  by 

1  -Yt(s)  0 

-Yt(s)  l+Yt2(S)  0 

0  0  1 

From  (11.3)  and  (11.4),  it  follows  that 

I  =  tr  Ct(t-s)  =  3  +  Yt(s)2, 


(11.3) 


[ Ct-1 ( t-s ) ]  = 


(11.4) 


and 


(11.5) 


II  *  tr  Ct-1 ( t-s )  =  3  +  Yt(s)2. 

Substituting  (H.5)  in  (8.4),  we  find  that  the  shear  stress  T<xy> 


is  given  by 


T<xy> 


/<- 

*  r\ 


G'  (s)  )</(Yt(s) )  {-Yt(s) }  ds. 


(11.6) 


where 


Uc  i_ 

^(Yt2s))  =  f2(3+Yt(s)2,  3+Y t ( s ) 2 )  -  fl(3+Yt(s)2,3+Yt(s)2) . 

d0  f 

Similarly  the  first  normal  stress  difference  N]_  =  T<yy>-T<xx> 

def 

and  the  second  normal  stress  difference  N2  =  t<xx>-T<zz>  are 


given  by 


Ni  -  I  (-G'(o))  \p  (Yt2(0))  Yt(o)2  do, 
J0 


(11.7) 


N2  =  /  (-G'(a))  fx(3+  Yt(a),  3+Yt2(a))  Y t2(a)  da. 

J0 

The  function  G  and  ^  in  (11.6)  and  (11.7)  are  chosen  by  Laun  [20] 
to  agree  with  experimental  data  obtained  from  measurements  of 
steady  viscosity  and  transient  response  to  step-strain.  In 
Narain  and  Joseph  [24],  it  was  shown  that  the  data  of  step-strain 
experiments  are  reliable  only  at  times  larger  than  a  charac¬ 
teristic  relaxation  time.  For  a  fluid  confined  between  two  flat 
parallel  plates,  one  of  the  plates  is  given  a  sudden  displace¬ 
ment,  and  if  the  effects  of  transients  die  quickly  then  the 
assumption  of  step-strain  gives  useful  asymptotic  results  [24]. 
Therefore  if 


,  for  t>0 
o 


0  for  t< 0 


(11.8) 


B 


r-v  ^yv')  Vv r>  ■r.v./.y. 


we  get 


Yt  ( s ) 


0  for  0<s<t 


-Y  for  s*t 
o 

daf-Y0[l-H(t-S) ] 


(11.9) 


Substitution  of  (11.8)  and  (11.9)  in  (11.6)  and  (11.7),  one 
arrives  at 


ij<xy> 


■/“ 


(  -G  *  (s)  )  tp  (Yn)  Y  0  ds. 


(11.10) 


and 


Ni 


•/' 


(  -g’(s)  )  <P  (Yq)  Yq2  ds  . 


(11.11) 


Now  (11.10)  yields  the  remarkable  result 

T<Xy>(t)/yo  =  ^(tJ/Yo2. 

Eqs.  (11.10)  and  (11.11)  are  in  agreement  with  the  semi -empirical 
formulas  (15)  and  (16)  of  Laun's  paper  [20].  The  modeling  for 
low  density  polyethylene  was  done  with  the  following  choices  of  G 
and  P  : 


G(s> 


=  Z  31  eXP  (  “  T_> 


(11.12) 


and 


P  ( Y 1 2  ( s )  )  =  fi  exp  (-n^  j  Yfc(s)  |  )  +  f2  exp  (-n2  I  Yt(s)  I  )- 
In  (11.12);  Ti,  i  =  1,2...  is  interpreted  as  a  discretization  of 


:--:y 

. 


•  *-  VJ 
>"v  w. 

m 


S». 


JRL. 


For  the  choice  of  constants  in  (11.12)  and  (11.13),  the  reader  is 
referred  to  Laun  [20].  The  remarkable  validity  of  eq.  (11.11) 
has  been  noted  by  Wagner  ([37],  [38]),  and  Laun  [20].  This  vali¬ 
dity  of  (11.11)  for  a  low  density  polyethylene  melt  is  shown  in 
Fig.  11.1. 

We  also  note  that  step  in  shear  rate,  start  up  of  elonga- 
tional  flows,  and  shear  dependent  steady  viscosities  for  melts 
can  all  be  modeled  (see  [19])  by  (8.4).  Successful  modelings 
done  by  Wagner  [36,  37],  Chang  and  Lodge  [7],  and  Laun  [20]  can 
be  reinterpreted  in  the  context  of  the  K-BKZ  model.  The 
superiority  of  this  model  in  various  flow  regimes  of  moderate 
strain  rates  has  already  been  emphasized  by  Tanner  [333. 


APPENDIX 


A)  Let  us  assume  the  kinematics  of  section  3  associated  with 

configuration  <  and  write 

S  =  V  (X,T  )  , 
fco 

where  (A.l) 

<(P)  =  X,  tQ  <  x  <  t. 

The  constitutive  theory  leading  to  (3.18)  assumes  that  the 
material  is  undistorted  up  to  time  tQ.  But  this  means  that  any 
other  undistorted  configuration  u  of  the  particle  P  at  a  time  t0 
could  at  most  be  a  rigid  rotation  of  the  configuration  K  .  To 
describe  this,  let 

K(P)  *  X, 

and  (A. 2) 

U(P)  =  Y. 

Let  X  be  the  invertible  function  which  maps  <  onto  U  such  that 

Y  =  A (X) .  (A. 3) 


Now,  for  t0  <  t  <  t,  we  introduce 


From  (A. 2) 


and  ( A . 3 ) : 


( *-i( Y)  ,t  ) 
(  Y  ,  I  )  . 


(A. 4) 


UF*  (t)  =  K F t  <  r  ) V  \  ,  (A. 5) 

~  '-o  ~  uo 

for  M  [tQ,  t  ]  .  In  (A. 4);  ‘Ft0< T  )  ,  Ft  (  T  )  and  V  *  are  respec¬ 

tive  gradients  of  aXt0  (  Y  ,  •■  )  ,  \t  (X,t),  and  x.  The  relation 
in  (A. 5)  implies  that  the  two  configurations  (not  pairs  of  par- 


tides)  K  and  u  undergo  the  same  changes  in  shape-history.  The 
two  configurations  <  and  u  are  undistorted  configurations  of  the 
same  material  element  only  if 


orth  (A. 6) 

is  an  orthogonal  tensor.  The  orthogonality  of  VA  in  (A. 6)  yields 
the  result 

K  li 

(A. 7) 


K  U 

yt<T)  -  ytcn 


for  all  tQ  <  x  <  t.  Now  the  natural  definition  of  the  response 


K»  I* 

functionals  K  and  ^  (as  introduced  in  (3.12))  are: 


<  X  #  K  K 

t  =  n  (  yt<T) ,  pto(t)> 


(A. 8) 


U  -v  #  u  Vi 

*  K  (  ot(n,  Mpto(t)). 

On  substituting  (A. 5)  and  (A. 7)  in  (A.8)^  we  find 

<  * 

T  -  V(  <  UUt<T>,  UFt  ft) (VA)"1) . 

~  '  ~  ° 

Since  ( V X ) -X  in  (a.9)  is  orthogonal,  use  of  the  definition  (3.14) 
of  isotropic  material  symmetry  of  the  configuration  <  implies 


(A.9) 


that 


T  -  K  Sri  (UUt(T  )  ,  MFt  (t)  (VA  )-1) 

'tfcV*  ~  -  ° 


(A. 10) 


= K  H  (uut<T) '  UFt0<t)) 


Therefore,  for  an  isotropic  simple  material,  combining  (A. 8)  and 
( A . 10 ) ,  we  have 


T  =  <  X  <MUt(  T)  ,  Ft  (t)  ) 

't^o 


(A. 11) 


-  u  ^rl  (uUt(T)  *  wFtft (t)) . 

The  result  above  establishes  that  for  an  isotropic  simple 
material ,  two  undistorted  real  states  have  the  same  represen¬ 
tation  for  their  response  functionals.  Now  if  the  arguments 
given  through  (3.12)  to  (3.18)  are  repeated,  we  find 

V 


where 


t  .'-y  [ct(t) ,  ct(tQ) ] 

-  -y  t ct ( t ) ,  ct(tQ)], 

k  y  def 

Ct ( T )  *  Ct(T)  *  Ct(T) 


(A. 12) 


for  tQ  <  t  <  t.  Furthermore j  if  the  material  is  undistorted  up 

Kjj 

to  time  tQ/  we  postulate  that  the  representation  of  rj  in  (3.18) 
must  be  such  that 

w.  *■ 

t  »  *y  ict(t) ,  ct(t0) i 

K.  i 

*  tct(T) ,  ct(t*n 

for  any  t*  <  tQ.  The  motivation  for  (A. 13)  is  the  simple 

K 

undistorted  rest  state  XtQ  (X,T)  =  x  for  all  t  <  tQ . 


(A. 13) 


B)  For  an  elastic  solid,  (3.5)  reduces  to 
T  ( x  ,t ) 


( B  .  1 ) 


*  Vftc(t))' 


and  therefore 


T(C,x)  =  4>1(Fto(T)  )  . 

In  (B.l),  *$1  such  that  K  (1)  =  0.  If  the  current 

configuration  is  used  with  the  help  of  (3.6),  we  can  define 


T  (  C  r T  >  =  $.(Pt(T)  Ft  ft)) 


(B.2) 


3.  $ 


2(Ft{x),  Fto(t)), 


whereK4>2!  sym.  Furthermore,  if  the  arguments  leading  to 

(3.18)  are  repeated  for  the  special  case  of  isotropic  elastic 
solid  given  by  (B.l)  and  (B.2),  it  is  easily  seen  that  there  is  a 


(B.3) 


response  function  <t>  such  that 

T<  ?,t)  =  K<t>  (Ct(T)  ,  Ct(t0)  )  , 

and  it  has  the  obvious  values 

T(X,t0)  *  K<f(Ct(t0)  ,  Ct(t0) ) 

=  0, 


T  ( x ,  t )  =  *  (1,  Ct(t0) ) . 

The  result  in  (A. 11)  of  Appendix  A  emphasizes  that  the  choice  of 

any  other  undistorted  configuration  will  yield  the  same  represen- 

£ 

tation  of  $  in  (B.3)  .  We  use  this  to  assert  that  there  is  a 
response  function  ^ct(t0)  (determined  by  the  natural  state  of 
the  elastic  solid)  such  that  the  Cauchy  stress  at  a  time  [  tQ ,t  ] 


is  given  by 


T(£,t>  =  $ct  ( t0 )  <Ctm).  (B  .4 ) 

Furthermore,  in  this  special  case,  the  Objectivity  requirement  in 


(3.20)  reduces  to 


Ct(tG)QT  (Q  Ct(T)QT>  =  Q  [*ct(t0>  (Ct(T))]QT 


(B  .  5 ) 


for  every  orthogonal  tensor  Q, 


KNOT 


,  .-‘.vV; 


7.yy, 


l  .v.v 


»  m  m 


HU 


.  /. 
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